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We construct generally applicable short-time perturbative expansions 
Qh. for some fidelities, such as the input-output fidelity, the entanglement fi- 



delity, and the average fidelity. Successive terms of these expansions yield 



characteristic times for the damping of the fidelities involving successive 



powers of the Hamiltonian. The second-order results, which represent the 
damping rates of the fidelities, are extensively discussed. As an interest- 
ing application of these expansions, we use them to study the spatially- 
correlated dissipation of quantum bits. Spatial correlations in the dissipa- 
tion are described by a correlation function. Explicit conditions are derived 
for independent decoherence and for collective decoherence. 
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In various fields of physics, the study of open quantum systems plays an 
important role. For example, decoherence was recognized as a major problem 
in realizing quantum computation [1,2]. Decoherence in quantum comp úters 
mainly results from the coupling of the quantum bit (qubits) to the environment. 
In general, It is not practical to look for exact solutions of these complicated 
systems, which consists of many spatially-correlated qubits interacting with a 
reservoir. 

In this paper, we propose a perturbative approach to the study of open quan- 
tum systems. We construct generally applicable short-time perturbative expan- 
sions for some fidelities, such as the input-output fidelity, the entanglement fi- 
delity, and the average fidelity. A similar perturbative expansion for coherence 
loss has been proposed in a recent paper [3], where the coherence loss is measured 
by the quantity 5 (t) = tr [p (í) — p 2 (í)], in which p (t) indicates the reduced den- 
sity of the system. Since 5 (t) has no direct physical meanings, in this paper, we 
choose the fidelities rather than 5 (í) as measures of decoherence. The fidelities 
are important quantities and have been widely used in quantum coding theory 
[4-7]. Through the perturbative expansions, we can get some general relations 
between the fidelities, such as the demonstration that the entanglement fidelity 
decays more rapidly than the average fidelity. 

As an interesting application of this perturbative approach, we use it to study 
spatially-correlated dissipation of the qubits. The latter question is of great 
practical importance in designing quantum computers. Spatial correlations in 
pure dephasing of the qubits have been analyzed in Refs. [8] and [9]. However, the 
calculations of spatial correlations in general dissipation of the qubits, including 
the dephasing and the relaxation, are far more involved. This question is hard to 
solve either by any existing exact approaches or by solving the master equation 
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[10]. Fortunately, here we show, through the perturbative expansions, spatial 
correlations in general dissipation of the qubits can be easily analyzed. Spatial 
correlations are described by a correlation function. From this result, we derive 
explicit conditions for independent decoherence and for collective decoherence, 
which are two important ideal circumstances in dissipation of the qubits. 

1 Perturbative expansions for the fidelities 

We consider an open quantum system, the total Hamiltonian of which is expressed 
as 

H T = H + Hj + H env , (1) 

where H and H env indicate the free Hamiltonians of the system and of the 
environment, respectively. Hj is the interaction Hamiltonian between the system 
and the environment. First, we suppose that the system is initially in a pure 
state l^o)- If there is no coupling between the system and the environment, at 
time t the system is in the state 

\y(t))=e- lHot / h \y ) = Uo(t)\y ). (2) 

But in reality, coupling of the system to the environment is inevitable. So at time 
t the system is in fact described by the reduced density 

P (t) = tr env [exp (-iH T t/h) p env (0) <g> |tf > (^o| exp (iH T t/k)} , (3) 

where p env (0) is the initial density of the environment. Decoherence of the system 
due to this inevitable coupling can therefore be measured by the fidelity between 
the state (2) and (3), which has the form 

F(t) = (V \U+(t)p(t)U (t)\y ). (4) 
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The input-output fidelity is first defined in Ref. [4]. Here the definition is slightly 
modified to Eq. (4) to make it more suitable for measuring decoherence of the 
system. 

In general, the reduced density p (t) in Eq. (3) is hard to calculate. So there 
is difficulty to obtain the fidelity (4). Fortunately, in practice, the short-time 
behaviors of the system under dissipation are of most interests. Hence we do not 
need to exactly calculate the fidelity (4), but expand it into a short-time power 
series 

F(í) = l-A_^_.... (5 ) 

Following Eqs. (2), (3) and (4), it is not difficult to obtain the expansion coefli- 
cients. Up to order t 2 they read explicitly 

1 = (6) 

|jr = (tfol tr env {[Hj, [Hj,p env (0) <g> |^ ) (tf |]] 

+ m, H },p env (0) ® |* > (* |]} |^o) (7) 

= ( H D S ,env-{( H l) 2 s) env i 

where the symbol (• • -) s env stands for the average value over the system and the 
environment, i.e., (H I ) S = (tf | #/ l^o) , and {Hj}^ env = tr env [p env (0) (ií/)J. It 
is obvious that the expansion coeflicients in Eq. (5) furnish characteristic times r n 
associated with decoherence processes involving H I to order n. The second-order 
coefficient \ is of special interest. It is the coefncient of the first non-trivial term 
in the expansion. Under the short-time approximation, ^ measures the damping 
rate of the fidelity. In most cases, it is sufficient to consider the expansion up to 
order t 2 . 

In the above, we have assumed that the system is initially in a pure state. 
The input-output fidelity (4) does not apply for the mixed states. For a mixed 



input state of the system, there are some fidelities defined, such as the entangle- 
ment fidelity and the average fidelity [5]. These two fidelities are widely used in 
quantum coding theory. Suppose \^ rs ) is a purification of the density p s , i.e., the 
state |\EVs) satisfies tr r (\^f rs ) fèrsl) — Ps, where the symbol r denotes an ancillary 
system. Similar to Eq. (4), the entanglement fidelity F e is defined as 
F e (t) = (tt„| C/ + (t) tr env [exp (-iH T t/k) p env (0) 

<g> |*„) (*„| exp (iH T t/h)\ U (t) |*„) (8) 




where F e (t) is furthermore subjected to a short-time power series expansion. 
In Ref. [5], it has been proven that the entanglement fidelity defined above is 
an intrinsic quantity of the system s, i.e., it does not depend on the specific 
purification \^ rs ). From Eq. (5), we easily obtain 

— = (9) 

¥t = (hi) -(mi) . (io) 

2T2 e ^ ' s,env \ ' s / env 

Comparing Eqs. (6) (7) with Eqs. (9) (10), we see, the characteristic times 
for the entanglement fidelity have the same forms as those for the input-output 
fidelity. Only change is that the symbol {Hï) s now means (Hj) s = tr s (p s H I ). 
So the entanglement fidelity is a natural extension of the input-output fidelity to 
include the mixed input states. Eqs. (9) and (10) also show, the characteristic 
times for the entanglement fidelity are all intrinsic properties of the system s. 
The initial density of the system can be expressed as a mixture of pure states, 

i.e., 

p a =Ettl*i><*il> ( n ) 
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where pi satisfy X) pi — 1. The average fidelity F a is defined as 

i 

F a (t) =E PiF 



(12) 



_ 2 t f 



where F denotes the input-output fidelity for the pure state Since the 

expression (11) for the density p s is not unique, unlike the entanglement fidelity, 
the average fidelity is not solely defined for a definite p s . From Eqs. (12) and 
(4), we get the characteristic times for the average fidelity 

— = 0, (13) 

Tia 

5" = -(?*<*'>?) ■ < 14 > 

ía \ » / env 

where the symbol (H I ) i indicates Hi Comparing Eq. (10) with Eq. 

(14), we get an interesting inequality. For any interaction Hamiltonians Hj, 
there is the operator inequality 

EMHi)ï> (EMHi)^ =(Hj)1 (15) 
Hence Eqs. (10 ) and (14) yield 

- > — ' ( 16 ) 

which suggests, the entanglement fidelity decays more rapidly than any average 
fidelities. In Ref. [5] it has been proven that the entanglement fidelity is always 
less than the average fidelity. Here we show this fact from another aspect. 

2 Spatially-correlated dissipation of the qubits 

As an interesting application of the perturbative approach to open quantum sys- 
tems, we consider a practical question: decoherence in quantum computers. This 
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decoherence is due to the inevitable coupling of the qubits to the external envi- 
ronment. A few papers [1,2,8,9,11-14] have been published on this subject with 
some simplifications, such as omitting spatial correlations in the decoherence [11- 
14], or omitting relaxation of the qubits [8,9], or omitting both of them [1,2]. 
However, in real circumstances, such as in the ion-trapped quantum computers 
[15], relaxation of the qubits has notable contributions to decoherence [11-13]. 
On the other hand, spatial correlation properties of decoherence play an impor- 
tant role in the choice of the decoherence-reducing strategies. For independent 
decoherence and for collective decoherence, the decoherence-reducing strategies 
are quite different. So here we consider a more practical model of decoherence. 
This decoherence is described by a spatially-correlated amplitude damping, which 
includes both the dephasing and the relaxation of the qubits. The environment 
is modelled by a bath of oscillators with infinite degrees of freedom and the mode 
functions of the bath field are chosen as plane waves. The Hamiltonian describing 
this decoherence process has the form 

H = H + K ( e ~ lkr ' a k + elkna t) ^ + W)] + E fakata k ) , (17) 

1=1 k k 

where the Pauli operator o\ represents the / qubit and a k is the annihilation 
operator of the both mode k. The symbol r; denotes the site of the / qubit, and 
Ai, A 2 and g k are coupling constants. H in Eq. (17) describes the free evolution 
and the infernal interaction of the qubits. The qubits, whether in memory or in 
quantum gate operations, can all be described by the Hamiltonian (17). 

The Hamiltonian (17) is very complicated and it is hard to find its exact 
solutions. Fortunately, with the perturbative approach developed in the previous 
section, this complex system can be easily treated. To analyze the decoherence, 
we use the perturbative expansion for the entanglement fidelity, which returns to 
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the input-output fidelity if the initial state of the qubits is pure. The environment 
is supposed in thermal equilibrium, i.e., the initial density p env (0) of the bath 
has the following form in the coherent representat ion 

Penv (0) =Q f d 2 a k 1 exp ( j \ ak ) ( afe | ? (ig) 



where the mean photon (or phonon) number 



(Nu*) = V 



Í^UJ k \ 



(19) 



With this density, and substituting the Hamiltonian (17) into Eq. (10), we obtain 
the decoherence rate — for this system 



4= E ^K-rfaXAA.A^,)., (20) 

2e h,/2=l 

where the operator 

^ = A 1( j£ + A 2 <j£ (i = 1,2) (21) 
and the spatial correlation function 



tt 2 (r h - r h ) = 2 E ||fffc| 2 cos [fc (r íx - r Í2 )]coth ^ 7 



(22) 



2k B T ) J 

Spatial correlation properties of the decoherence is completely determined by the 
correlation function fl 2 (r^ — rj 2 ). 

In the following, we discuss two important ideal circumstances for the corre- 
lation function. Eq. (22) can be rewritten as 

Ü 2 (r h - ri2 ) = xJ2f(k) cos [k (r h - r h )\ , (23) 
k 

where / (k) is a normalized distribution satisfying J2 f (k) = 1 and x is the 

fe 

normalization constant 



= 2 J2 M 2 coth ^ 



x = 2Ek! Z coth(^V (24) 



The expression of / (k) is given by comparing (23) with (22). Its explícit form 
depends on the coupling coeficient \gk\ , whereas the latter is determined by the 
specific physical model for quantum computers. As a simplification, here we as- 
sume that the distribution / (k) can be approximated by a Gaussian function with 
an expectation value k and a variance Ak, respectively. With this simplification, 
Eq. (23) reduces to 



exp 



-\{Akf{r h -r l2 f 



(25) 



fi (r h -r h ) w x cos \ k (r h - r h ) 
Suppose d is distance between the adjacent qubits. If d satisfies 

(Ajfe)d»l, (26) 

Eq. (25) yields 

(r h ~ r l2 ) ^ x5 hh , (27) 
and from Eq. (20) the decoherence rate is simplified to 

± =£ x ({AAtf) (28) 

T 2e 1=1 

which suggests, the total decoherence rate of L qubits equals the sum of the 
decoherence rates of individual qubits. So in this circumstance, the qubits are 
decohered independently. Eq. (26) is the condition for independent decoherence. 
Most of the existing quantum error correction schemes are designed to correct 
for the errors induced by independent decoherence [16-24]. 

Apart from the independent decoherence, there is another ideal circumstance 
for the spatial correlation function. Suppose there are 2L qubits. Two adjacent 
qubits make up a qubit-pair. So we have L qubit-pairs. The two qubits in the 
l qubit-pair are indicated by l and t , respectively. If distance d between the 
adjacent qubits satisfies the condition 

kd « 1 and (Ak) d « 1, (29) 



form Eq. (25) the spatial correlation function remains a constant for the two 
qubits in each qubit-pair. The decoherence rate ^- for L qubit-pairs thus becomes 

4 = E n 2 (r h - r h ) (A (A h + A f ) A + AA) . (30) 

T 2e I 1 ,ia=l 

So under the condition (29) the two qubits in each qubit-pair are decohered col- 
lectively. In the collective decoherence, the decoherence rate is sensitive to the 
type of the initial states. If the initial state of the qubit-pairs is a co-eigenstate of 
all the operators Ai + A[, form Eq. (30), the second-order decoherence rate ^- re- 
duces to zero. With these states decoherence of the qubits can therefore be much 
reduced. The co-eigenstates of the operators A\ + A\ are called the subdecoherent 
states. In fact, an arbitrary input state of L qubits can be transformed into the 
corresponding subdecoherent state of L qubit-pairs by the following encoding 

1-1) - 1-1, +1), (31) 
1 + 1) -1+1,-1), { } 

where | — 1) and | + 1) are two eigenstates of the operator A\. Obviously, the 
encoded state is a co-eigenstate of the operators A { + A\ with the eigenvalue 
0. This encoding has been mentioned in [8] and [25] and extended in [26] and 
[27] to reduce decoherence in general circumstances. Here we derive the working 
condition (29) for this encoding. 

In the above, k and A/c are introduced phenomenologically. It is assumed that 
we have no knowledge about the coupling coemcient \gk\ ■ There is an interesting 
case in which the spatial correlation function can be calculated exactly. Consider 
one dimensional quantum comp úters. In the continuum limit, it can be assumed 
that 

£ M'···oc / dw k u k e-^ ■ ■ ■ (32) 

where u c is the cut-off frequency whose specific value depends on the particular 
nature of the physical qubit under investigation. The form (32) of \gk\ 2 was also 
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used in Refs. [1] and [8]. We consider two circumstances. In the high temperature 
limit, i.e., T >> the correlation function (22) is simplified to 

íl 2 (r h - r h ) oc J- du k e-»*/«° cos [<f (r h - r h ) 

(33) 

= Ak B T uj 2 c 

^>c í+^r^-r^ )/»]*' 

where v indicates velocity of the noise field. If distance d between the adjacent 
qubits satisfies d » v/u c , Eq. (33) tends to a delta-function and the qubits are 
therefore decohered independently. On the other hand, if d « v/cu c , we have 
fl 2 (d) = Q 2 (0) and the adjacent qubits are then decohered collectively. In the 
low temperature limit, things are much similar. If T << coth (^^) ~ 1 
and Eq. (22) yields 

n 2 (r h - r l2 ) oc 2 / °° du k u k e~^ cos (r h - r h ) 



l+w§(ri 1 ~r l2 f /v 2 

Eq. (34) suggests, the above conditions for independent decoherence and for 
collective decoherence still hold at low temperature. The type of decoherence is 
mainly determined by distance between the adjacent qubits and by the cut-off 
frequency. Temperature of the environment hardly influences the decoherence 
type, though it determines the decoherence rate. 

3 Summary 

In this paper, we develop short-time perturbative expansions for some widely- 
used fidelities. From the expansions, we demonstrate some interesting relations 
between the fidelities. Perturbative expansions for the fidelities can be used for 
studying open quantum systems. As an example, we consider spatially-correlated 
dissipation of the qubits in the quantum computer. Spatial correlations in the 
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dissipation are described by a correlation function and we successfully derive 
the explícit conditions for independent decoherence and for collective decoher- 
ence. This example suggests, the perturbative expansions for the fidelities may 
be proven as a useful tool for studying open quantum systems. 
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